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We investigate effects of nuclear medium on antiquark distribution in nuclei applying the results 
of a recently developed relativistically-covariant self-consistent model for the pion and the isobar. 
We take into account Fermi motion including Pauli blocking and binding effects on the nucleons 
and medium effects on the isobar and pion leading to modest enhancement of the pion light-cone- 
momentum distribution in large nuclei. As a consequence the Drell-Yan cross-section ratio with 
respect to the deuteron exceeds one only for small values of the light-cone momentum. 



I. INTRODUCTION 

Improving our understanding of the quark and gluon 
degrees of freedom of nucleons bound in nuclei necessi- 
tates further efforts in spite of numerous investigations 
and successes The planned new nuclear Drell-Yan 
(DY) scattering experiment SeaQuest at Fermilab (E- 
906) strongly motivates recent advancements in 

analysis of nuclear parton distributions @ as well as 
attempts directed at more reliable and accurate calcu- 
lation of experimentally observable quantities. In line 
with the second objective is our aim to update and ex- 
tend to larger x values previous calculations of nuclear 
DY effect presented in Refs. @, 0. Furthermore, the 
physics interest is twofold: (i) to investigate the role of 
non-perturbative physics (pion cloud) in the u — d asym- 
metry in the nucleon, and (ii) to study possible anti- 
quark enhancement due to presence of virtual mesons in 
the nuclear medium. From symmetry properties of QCD 
we may infer that pion is one possible source of non- 
perturbative quark-antiquark physics, in particular u — d 
asymmetry. 

In conventional models based upon meson exchange 
nuclear binding of nuclear matter comes for about 50% 
from (virtual) pions present in the nucleus. Can we ob- 
serve these? Some indication for pions is present in EMC 
effect enhancement around x = 0.1, which can be as- 
cribed to the fact that pions (and heavier mesons) carry 
a fraction of the momentum sum rule. Can one see these 
pions more explicitly, e.g. in the form of an enhancement 
of anti-quarks in the nucleus? Anti-quarks can be probed 
directly in Drell-Yan scattering but previous experiments 
Q within the experimental uncertainty of about 10% did 
not show a nuclear enhancement; results of calculations 
varied strongly. 

In practice one can distinguish two main types of the- 
oretical interpretations of the classical EMC effect: (i) 
in terms of nucleon as constituents which are bound but 
not modified in the medium, and (ii) in terms of off-shell 
nucleons with medium modified structure functions, e.g. 
through scalar and vector fields acting on the quarks. 
In the first category one has the non-relativistic models 



which use a computed spectral function that accounts 
for large removal energy (50 MeV in nuclear matter) and 
Fermi motion due to correlations. This approach can re- 
produce the observed slope of the reduction for x < 0.5 
in the EMC ratio 2F^{x) / AF$(x) but not the behavior 
around the minimum in the ratio at x — 0.8. The latter 
seems to require ad-hoc off-shell effects @. In the sec- 
ond category [H|[llj]) one usually starts from the Walecka 
model in the mean-field approximation which has a small 
net binding effect (8 MeV per nucleon) and hence yields a 
very small EMC effect, and then one adds the effect of ex- 
ternal scalar and vector fields. Since in the present study 
we are interested in the antiquark distribution for x < 0.4 
we rely on the conventional convolution approach using 
a parameterized nucleon distribution, which has the two 
above mentioned parameters related to the removal en- 
ergy r\ and the Fermi momentum. 

In Sec. II we study the antiquark distributions in the 
free nucleon and determine the off-shell 7rNN and 7rNA 
form-factors which lead to good description of the isovec- 
tor part of the proton antiquark distribution by the pion 
cloud. In Sec. Ill we turn to general discussion of the 
nuclear effects consisting of binding and Fermi motion 
of nucleons and modification of the pion cloud. Detailed 
consideration of the medium effects on the nucleon's pion 
cloud follows in Sec. IV where expressions are derived for 
the pion light-cone distributions originating from the 7rN 
and 7rA states. We emphasize the careful treatment of 
the in-medium delta baryon based on a complete rela- 
tivistically covariant basis for its dressed propagator. Nu- 
merical results for the in-medium pion distribution and 
DY cross-section ratios for nuclear targets relative to the 
deuteron are presented and discussed in Sec. V. Finally, 
Sec. VI contains a summary of our results. 



II. ANTIQUARKS IN FREE NUCLEONS 

Before turning to the nuclear case we want to inves- 
tigate whether the pion cloud approach we will use in 
the medium can reproduce the observed flavor asym- 
metry in the free nucleon. The distribution of anti- 
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quarks in the nucleon can be decomposed into a flavor- 
symmetric isoscalar part (originating from gluon splitting 
and possibly meson cloud) and a non-perturbative isovec- 
tor meson-cloud contribution. The latter can be consid- 
ered to be the source of the u — d asymmetry. In addi- 
tion to the pion cloud of the nucleon we include also the 
isobar with its pion cloud since it was shown to give sub- 
stantial contributions [12j . In this way one can constrain 
the 7rNN and 7rNA form factors using the empirical an- 
tiquark u-d flavor asymmetry. The physical free nucleon 
state is expressed approximately as 



(1) 



\N) = VZ\N) baxe + a\Nir) + p\Air). 



Neglecting off-shell effects the light-cone momentum dis- 
tribution of a quark with flavor / in a proton can be 
written as (B = N, A): 

q f (x) = Zqf, haxe (x) + c * [ -r-f Bz/N (y)<if,b^ c ( x /y) 
d ■ lj x y 



+ I -r /N (y)q^(x/y) 
y 



(2) 



where Ci (i labels the charge states) are the appropriate 
isospin Clebsch-Gordan coefficients, qf' harc (x) is the par- 
ton distribution in the bare Bi baryon and q^(x) is the 
pion parton distribution function. 

Attributing the asymmetry in the u and d antiquark 
distributions to the nucleon meson cloud we are con- 
cerned with the pion light-cone distribution in the nu- 
cleon which gets contributions from final states with ei- 
ther nucleon or isobar: 



r /N (y) = r N/N (y) + r A/N (.y). 

The nucleon term was calculated by Sullivan 13]: 
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with y = (fco + k^)/M being the pion light-cone momen- 
tum fraction, with M the physical mass of the nucleon (as 

a convenient scale), F^ N (t) is the off-shell form- factor 
of the ttNN vertex, while g^NN is the tt°NN coupling. 
The free-pion propagator, D®, appears in the above ex- 
pression in the form (t + m 2 )" 1 , where t = —q 2 where 
q denotes the pion four-momentum. The isobar contri- 



bution also plays an important role [12j despite the kine- 
matical suppression coming from the isobar-nucleon mass 
difference. In Ref. [12] it was calculated using the free 
isobar propagator, i.e. neglecting its width. A complete 
rclativistically covariant treatment of the isobar in vac- 
uum and nuclear medium was introduced in Ref. [l4 | and 
we use that formalism to take into account the vacuum 
width consistent with the measured pion-nucleon scat- 
tering phase shift in the spin-3/2 isospin-3/2 channel. 
The full Lorentz structure of the vacuum propagator of 
the Rarita-Schwinger field can be expressed in terms of 



10 Lorentz scalar functions II 4, 1 1 511 which contains both 
spin-3/2 and spin-1/2 sectors [lg. However, using the 
convenient basis from Ref. [l4| it turns out that a single 
term, namely the (on mass shell) positive energy spin- 
3/2 contribution gives the dominant contribution and all 
others (some terms in the propagator are identically zero) 
are completely negligible. In the notation of Ref. [lj] this 
is the coefficient of the projector sum Qffu = Q[i%] + -^[55] 

which we denote by G^(p). The pion light-cone distri- 
bution originating from the Att state then can be ex- 
pressed as: 
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where p and p ' are the four-momenta of the nucleon and 
isobar, t = -(p- F , ^ . rrNA 
the 7r7VA vertex. The form factor 



p') 2 and are the form factors of 



p 2 - (M + m^f 



A 2 



(6) 



with A = 0.97 GeV (and g^NA = 20.2 GeV r ) was used 
in Ref. [l5[ and shown to give a good fit to the relevant 
pion-nucleon phase shift. For the irNN and irNA off- 
shell form factors which take into account the off-shell 
pion we take a dipole form: 
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with X standing for N or A. 

In order to calculate the d — u distribution for the free 
proton we assume that the pion sea is isospin symmet- 
ric leaving only the contribution of valence distributions. 
The final state with nucleon contributes through the pres- 
ence of 7r + with distribution 2f n N / N y while the isobar 
final state can have also a ir + with isospin weight 1/3 
or a 7r~ with minus sign (because of the pion valence u 
distribution) with respect to f* A / N giving in total: 



(d-u) p (x) 



y 



{y)-\r' A/N {y))ql{x/y), 



(8) 

with q*(x) denoting the valence parton distribution of 
charged pion. In Fig. [T] we show the pion distributions 
with nucleon final state (solid line) and isobar final state 
(dash line). For the form- factor cutoff we used the follow- 
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0.95 GeV and A^ 



0.75 GeV. The 



ing values: A;^ 
bare-nucleon probability then takes the value Z = 0.69 
which suggests that higher-order terms with more than 
one pion do not contribute significantly. The calculated 
value for the d—u asymmetry in the free proton is shown 
in Fig. [5] by solid line and compared to the result using 
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FIG. 1. (Color online) Pion distribution in the free proton: 
r° N/N {y) shown by solid line and r~ A/N (y) by dash line. 



the u and d fits CT10 [17] (dot line). Also shown are 
separately the contributions from the nucleon final state 



(dash line) and isobar final state (dash-dot line). These 
results are quite similar to the d— u obtained in Ref. [l8[ , 
although there the infinite-momentum-frame formalism 
was used with suitably adjusted values of the 7rN and 
7rA form factors. 



III. NUCLEAR EFFECTS 

First calculations of the nuclear Drell-Yan process 
suggested an enhancement coming from the medium 
modification of the pion cloud. However, the experi- 
mental data [§| did not show that enhancement within 
a 10% uncertainty. Later on other groups reported more 
detailed calculations of the Drell-Yan ratio with a large 
variation in results as shown in Refs. [MM- Here we con- 
sider the ratio of the cross sections of proton- nucleus and 
proton-deuteron scattering, 
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where A denotes the nucleus and its nucleon number. We 
specialize for the case of isoscalar targets for which the 
cross-section ratio becomes: 
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In the case when x\ is large, say x\ > 0.3 the second term 
in the numerator becomes negligible and only medium 
effect on the antiquarks plays a role. 

The (anti-)quark distribution in the medium can be 
modified in two ways, (i) through Fermi motion and bind- 
ing of the nucleon, and (ii) modification of the nucleon's 
pion cloud. To establish the connection to the (anti- 
)quark distribution of the free nucleon we use Eq. ([5]) 
which for the free proton gives (with isobar terms not 
written out for brevity): 

q p f (x) = Zf fMre (x) 1 ' ' 



l f N/N {y) 
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• 2«y (x/m 



Z = l- dyf N / N (y) = l-3 dyr° N / N (y), (12) 
Jo Jo 

where the last equality expresses flavor-charge conserva- 
tion. Similarly, the quark distribution for the nuclear 
proton can be written: 

q p (x) = Z A ql haTe {x) + \ [ A d -±f N l A {y) k barc (*/y) + 

° Jx y L 



2?/,barcO/2/)] + 



^r° N/A (y) Uf{x/y) + 2qf(xm}) 
y L i 



where isospin-symmetric nuclear medium was assumed. 
Adding the difference of the left-hand side and right-hand 
side of (fTTj) to the right-hand side of (fTB"]) and repeating 
the same procedure for the neutron we obtain 

? f (x) + q](x) = q p f (x) + q p fMie (x) f f N ' N (y)dy - 



1 dy 



f N/N (y)i p f , b!a My)-<ffMreW / f N/A (y)dy + 



dy 



f N/A (y)q p fM re(x/y) + (p^">) 



,; A djL 

y 



r° N/A (y)-r° N/N (y) 



■ q lo(x/y) + q l+(x/y) + q l-(x/y) 

where we used that 



(14) 



Z A = 1- 



dyf N/A (y). 



To proceed with the above expression one needs the bare 
antiquark distributions which could be determined from 
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FIG. 2. (Color online) The pion-cloud result for the d — u 
asymmetry in free proton (solid line) compared to the differ- 
ence of the proton d and u distributions (dot line) from the fit 
CT10 and data points by Fermilab E866/NuSea Collabo- 
ration [20]. The isobar contribution (dash-dot line) is negative 
and much smaller than the nucleon term (dash line). Used 



parameter values: = 0.95 GeV and A^ = 0.75 GeV. 



the antiquark distribution of the free "isoscalar" nucleon 
for two used pion parameter sets. For the distribution of 
the bare nucleon in the free nucleon we use the relation 
f N / N [y) — 3/ 7r N ' N (l — y) which follows from the proba- 
bilistic interpretation of these functions [18j and its ana- 
log for the in-medium case, f N ^ A {y) = 3/" N I A (1 — y). 
We remark that the latter relationship can only be ap- 
proximate since the support of the pion in- medium distri- 
bution is not strictly limited by value one, but in view of 
the similarity of the pion distributions in the two cases it 
should be a reasonable approximation for the estimate of 
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Eq. (jllj) and its analog for the neutron. A much sim- 
pler, though approximate procedure is just to subtract 
the meson-cloud contribution form the antiquark distri- 
bution of the physical nucleon. Indeed, using the fact 
that antiquark distributions at small x behave as X/x 
one can confirm that the first and second term on the 
right-hand side of (fTT|) combine to give the bare distri- 
bution if one adds to (fTTj) its neutron analog. Using the 
same argument about the small x behavior of antiquark 
distributions one can establish an approximate cancella- 
tion of the second and third as well as the fourth and 
fifth terms on the right-hand side of Eq. (fT4|) and of the 
corresponding terms involving the neutron. This simpli- 
fication was used in our previous work [6(, but in the 
present calculation we want to take into account these 
contributions with the bare antiquark distributions de- 
termined by the above mentioned subtraction of the pion 
contribution from the physical distribution. In Fig. [3] we 
show the proton-neutron average of the sum of the sec- 
ond, third, fourth and fifth terms in Eq. (TT4")) divided by 



FIG. 3. (Color online) Sum of the second, third, fourth and 
fifth terms in Eq. (|14p devided by the antiquark distribution 
of the free "isoscalar" nucleon. Solid line is for the pion pa- 
rameter set (1) and dash line for parameter set (2). 

a small effect. We observe that this contribution is indeed 
quite small as expected from the form of the nucleon an- 
tiquark distribution and pion (as well as related nucleon) 
light-cone-momentum distributions. The parameter set 
(1) is: M* = 0.89 GeV, E^ = 0, E^ = -0.1 GeV, E^ = 
0, g' u = 0.9, g' 12 = 0.3, 522 = 0.3, while the set (2) is given 
by: M* = 0.89 GeV, = 0, = -0.05 GeV, = 
Q,g'n — 1.0, g' 12 — 0.4, g' 22 — 0.3; where M* is the mean- 
field shifted nucleon mass, Ejy- the energy shift of the 
nucleon, and the delta's mean-field shifts and g'^ 
the Migdal four-fermion interaction parameters. 

Introducing a shorter notation for the bare nucleon 
contribution to the in-medium antiquark distribution: 



y 

<ff M r^)i f N/A (y)dy+ I ff N/A (yW fMTe (x/y) + (p^n), (15) 
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we can finally write the sum of in-medium proton and neutron antiquark distribution as 

<f/ A (x) + q f A (x) = £ &f&( 3 ,)$g n {x/v) + [r° /A (v) - r" /N (y)} [iU-/y) + 



oi-(x/y) 



(16) 
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FIG. 4. (Color online) The ratio of Fi(x) per nucleon for 
isospin symmetric nuclear medium with parameters: Fermi 
momentum pf = 250 MeV, r\ = 0.97, and deuteron. Pion 
contribution with parameter sets (1) and (2) is included in 
the results shown by solid and dash-dot lines. 



The convolution with fp b (z) takes into account Fermi 
motion and binding effects on the in-medium nucleons. 
For the function fp b (z) we take the result of Birse (2lj : 



V\)- 



(17) 



where e = pf/M, rj is a parameter with value slightly 
below one which takes into account the nuclear binding 
and Q(x) is the unit step function. In Fig.2]we show the 
ratio of the F% (x) structure functions for the isospin sym- 
metric nuclear matter and the deuteron. We assume neg- 
ligible medium effects in the deuteron and for the nuclear 
medium use the convolution model with light-cone distri- 
bution p7|) and parameters pf = 250 MeV and r\ = 0.97. 
In this way one can reproduce the negative slope in the 
classical EMC effect for 0.1 < x < 0.5 as shown for exam- 
ple in Ref . [22j ■ The experimental enhancement observed 
around x = 0.1 can be attributed to the pion enhance- 
ment as shown in the figure for the two parameter sets 
(1) and (2) used also for the plots in Fig. [3] and given 
above. The pion enhancement term was calculated by 
the convolution of the in-medium pion light-cone distri- 
bution enhancement relative to the free nucleon and the 
pion i*2 distribution. Note that expected shadowing ef- 
fects would lead to decrease of the nuclear cross section 
for x < 0.05. 








FIG. 5. Types of diagrams contributing to the in-medium 
pion distribution with outgoing nucleon. The same types ap- 
pear also with outgoing delta baryon. The dash line denotes 
the dressed pion propagator, while the solid and the double 
line correspond to nucleon and delta. 



IV. ANTIQUARKS IN BOUND NUCLEONS: 
PION CONTRIBUTION 



We now turn to consideration of the pion contribution 
to antiquark distributions in nucleons bound in large nu- 
clei which we model by an infinite system with appropri- 
ate average nuclear density. 

For corresponding pion properties in the nuclear mat- 
ter we use the results of a fully covariant self-consistent 
model developed in Ref. [23J. Compared to the case of 
the free nucleon the nuclear environment changes the 
pion propagator appearing in the Sullivan formula (£[]) 
and rcnormalizes the 7rNN as well as the 7rNA ver- 
tices through nucleon-nucleon correlations modeled by 
the Migdal four-fermion interactions [24| • Nucleon prop- 
erties are also affected and we take into account the bind- 
ing effects through mean-field mass and energy shifts con- 
sistent with approach in Ref. [23j. 

The inclusion of the dressed pion propagator is 
straightforward but the dressing of the 7rNN and the 
7rNA vertices requires summation of nucleon-hole and 
delta-hole bubbles. The types of relevant diagrams are 
shown in Fig. [5] For resummation of these diagrams we 
use the rclativistically covariant formalism introduced in 
Ref. [Ill and applied for pion and delta self energy cal- 
culation in Ref. [23[ . In the present case it concerns a 
different type of contribution which for the nucleon in 
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the final state can be written: 

K N = u(p') l5 ^u(p) ■ n^(g)g", (18) 

where q = p — p' is the pion 4-momentum, u(p) is the nu- 
cleon in-medium spinor (with mean-field shifts of mass 
and energy) and n M „ (q) is the resummed contribution of 
nucleon-hole and delta-hole loops. Using the decomposi- 
tion of nucleon-hole and delta- hole loops [23[ : 

n<S fc, (9)= E n^(g)Lfe)(g) + n£%) 1^9) 

where L^l (q) and T llll {q) have projector properties, 
I 



the nonzero contribution involving one particle-hole 
(nucleon-hole or delta-hole) loop comes from 

n« = sir (<%) 4" } (<z) + n^to 

where 5'n, 512,522 are the usual Migdal parameters [23[ 
where index 1 refers to the nucleon and index 2 to the 
delta and we took into account that: 

L$\q) ■ q» = L$\q) ■ q» = T^{q) ■ g» = 0. (21) 

In order to perform the summation involving arbitrary 
number of nucleon-hole or delta- hole loops it is conve- 
nient to introduce the following matrices [25[: 



fg'u g'i 2 

5n 

5l2 522 

\0 512 522 



/ n<? h) < h) \ 



5l2 



, n< L > = 



Xi 21 


n 














n (Ah) n (Aft) 

n (A? l ) n (A/») 
11 21 ix 22 



(22) 



The lowest order contribution ([2U)l can then be written as: 

ng = [(nW fl W) 11 + (nW fl W) 81 ]xa i) (fl) 
+ [(n( i )5 (L) )2i + (n^5 (L) )4i]4i 1) (g). 



(23) 



Higher order terms are accounted for by taking appropriate matrix elements of products of (H( L ^g( L ') matrices 
and the summation of terms with arbitrary number of loops is simply achieved by replacing (nWj^) in (|23p by 
n( L y L ) (1 - U^g^l)- 1 leading to: 



n 



(i _ nWflW)- 1 ]!! + [nw 9 w (i - nWjWj-V) £g, 1} (<0 

+ ([Tl^gW (1 _ nW^)^ + [Tl^gW (1 - iK^)- 1 ]^) Lg>( ff ). 



(24) 



Adding 5 M „ to n Mt/ in (|18[) we obtain the full contribution of the diagram with nucleon final state. Squaring its 
absolute value and performing summation over the spin projections of the nucleon in the final state and averaging for 
the nucleon in the initial state we obtain: 

~ \ R n\ 2 = A qq ■ (2(p - E» -qtf- £» -q-q 2 [M 2 + (p' - ■ (p - £»]) 

s,s' 

+ 2A qu ■ ({p - £» -nip'- S» -q+ip'- ■ u (p - E v N u) ■ q 

- q-u[M* + (p'-X v N u)-(p-Z v N u)}) 

+ A uu • (2(p - £» -nip'- £» • u - M 2 - (p' - £» ■ (p - £») , 



(25) 



E» r and £» r are 



where u is the 4- velocity of the medium (implicitly present also in Eqs. (|T8)) — (j2~4")) ), M* = Mjv + ^jv> ^jv allu ^jv 
the nucleon mean-field mass and energy shifts. The factors A qq , A qu , A uu are given by: 

2 



A qq — 2 



A n 



A — 9 



1 + A + q ■ u B/yJq 2 - (q ■ u) 2 

2 Re 



<fB/y/q* - (q ■ uf{\ + A~q -uB/y/q? - (q ■ U ) 2 ) 
2 



q 2 B/^qZ-{q-uf 



(26) 
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where the bar denotes complex conjugation, 

a ee (i - nw^)- 1 ^ + p^jW (i - nw^))- 1 ^, 

b ee [n( L y L > (i - n^pW)-^! + [n( L y L > (i - nW^))- 1 ]^ (27) 



and we used that: 
L$\q)q v 



q ■ u 



\Jq 2 -{q- u) 2 \/q 2 -{q- u) 2 

To compute the pion light-cone distribution per nu- 
cleon in the medium we integrate over incoming nucle- 
ons in the Fermi sea and outgoing ones above the Fermi 
sea, restricting the pion light-cone momentum fraction to 
the specified value by inserting a delta function and fi- 
nally divide by the nucleon density. The final expression 
obtained is: 



f N (y) = 3My 
p' x dp' ± 



fN 

2tt 



1 



32tt 3 P 3 f 
1 



dp3 







s,s' 



PF 

2 



ttNN\ 



where D 7r (q) is the in-medium dressed pion propa- 
g ator, b ee My - p 3 - ^M 2 +p 2 +pl, pf" n = 

yj 2by/M 2 +p 2 F - M 2 - b 2 , <& is the angle between pj_ 

and p'j_, and the irNN form factor F^j N (—q 2 ) was also 
included. 

The contribution coming from the delta baryon in the 
final state is made more involved by the complicated 
structure of the in-medium delta propagator |14j . How- 
ever, considerable simplification can be achieved by in- 
cluding only the two dominant contributions in the con- 
venient relativistically covariant decomposition since the 
imaginary part of the other components is typically two 
orders of magnitude smaller at nuclear saturation den- 
sity [26( | . The dominant contributions come from the 
P±dp± an d terms which were degenerate in the free 

delta case, but are different in the medium [141] . Sum- 
mation of particle-hole loops dressing the irNA vertex is 
analogous to the irNN case, with only difference being in 
the relevant matrix elements of li^g^ (l—H( L )g( L ))~ 1 , 
(29) replacing the expression ([2"4")) with: 



= (i - nW fl W)-% + (i - n<W]») ^\ q ) 

+ ([U^gW (1 - U^gWy'U + \n.MgM (1 - n^)" 1 ^) L$\q). 
The expression analogous to ([23]) in this case takes the form: 

I^aI 2 = \ Tr [{p - + MJImG^tf)] ( 9lia g v p + II* S*) gV, 



(30) 



(31) 



where Im G fiu (p') denotes the imaginary part of the in-medium A propagator for which we take the dominant contri- 
bution given in the basis used in Ref. [14| by just two terms: 



In this way the expression ([3"T]) takes the form: 

\ Y I^aI 2 = \A£> c») + 2A#> cg> + A® cW)l ImG^') 



(32) 



4(A) (P) , 04(A) JP) , 4 (A) (P) 



(33) 



where the expressions for A^" 1 , ^4 g „ , Ala , c g f , ci$) , ciu' 1 , , c^P , clu" 1 are given in the Appendix. The pion light- 
cone momentum distribution stemming from the process with the nucleon emitting a pion and a delta baryon is 
analogous to expression ([29]) and reads: 



f A (y) = MIy 



1 



TO^r J 327T 4 p^ j _ pp 
2tt 







d P3 



/•00 /•Z7T 1 

Jo Jo z a „, 



(34) 
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We checked by explicit numerical calculation that both 
Eqs. (|29|) and ([34]) have the correct low-density limit, i.e. 
reproduce the free nucleon and delta results. 



V. NUMERICAL RESULTS AND DISCUSSION 

For the computation of in-medium pion and isobar 
properties we rely on the recently developed relativisti- 
cally covariant self-consistent model presented in Ref. [23[ 
and used for the nuclear photoabsorption calculation in 
the isobar region in Ref. [271 ] . For the medium computa- 
tion we use the same 7rNN and 7rNA form factors as in the 
vacuum one and include them in the model of Ref. [23j . 
The values of the Migdal g' parameters which model the 
short-range nucleon and isobar correlations were taken 
in the range preferred by the results of Ref. [27J ; in this 
work a good description of the nuclear photo-absorption 
cross section in the isobar region was obtained. Binding 
effects for the nucleon are taken into account by the ef- 



0.30 

0.25 



- free nucleon 



M*=0.89GeV, g'„=0.8, g' 12 =0.3, g' 22 =0.3 

• - M*=0.85GeV, g',,-1 .0, g' 12 =0.4, g' 22 =0.3 
M*=0.89GeV, g', =1.0, g' 12 =0.4, g' 22 =0.3 




FIG. 6. (Color online) Nucleon contribution to the pion 
distribution (/" N ' (y)) for the in-medium nucleon compared 



fective (mean-field) mass M* and the energy shift £^. A to the case of free nucleon. 
consequence of the use of the mean-held approximation 
is a reduction of the in-medium pion distribution coming 
from the nucleon final state. Namely the dominant con- 
tribution to it is the term proportional to A qq in Eq. (|25[) 
with: 



(2(p - £» ■qip'- E» • q - q 2 [M 2 + (p' - £» 



■(p-£»])=-2M,y 



(35) 



which is the same expression as for the free nucleon, ex- 
cept that M* appears instead of M. 

Since M*/M < 1 a further suppression in addition to 
that from the Pauli blocking is obtained, depending on 
the actual value of M*/M. The latter is difficult to con- 
strain since observables generally are only sensitive to the 
combination M* + E]^. Since our aim is to make com- 
parison with experiments on finite nuclei (rather than 
nuclear matter) with an average density smaller than the 
saturation density we assume small values for the energy 
shift in the range zero to = 0.04GeV, corresponding 
to effective mass values in the range of of 0.85 GeV and 
0.89 GeV. These values are close to the ones used in more 
elaborate treatments of nuclear matter [28|,[29| where val- 
ues of 0.8 0.85 GeV at saturation density give good 

agreement with observables. The mean-field shifts of the 
isobar mass and energy are chosen in such a way so that 
they re pro duce the isobar-nucleon mass difference used 
in Ref. [27| . This means E^ = -0.05 GeV and -0.1 GeV 
and zero for the energy shift. 

In Figs. [5] and [7] we show the pion distributions 
N/A^ anc j jtt a/a^ f or in-medium nucleons for dif- 
ferent parameter sets. For the nucleonic distribution one 
observes a reduction coming partly from the Pauli block- 
ing of the nucleons in the medium and partly from the M* 
effect (which leads to a suppression roughly by the factor 
(M*/M) 2 ). The pion broadening in the nuclear medium 
only partly compensates these effects and a net reduction 



is the result. This is not completely surprising since the 
computations of Ref. 23j do not lead to appreciable soft- 
ening of the pion spectrum in the medium which would 
result in enhanced pion distribution. In this respect the 
pion dressing of Ref. 12311 i s not significantly different 
from an older calculation 3u| which used a nonrelativistic 
treatment of the isobar and a softer pion-nucleon-delta 
form factor. On the other hand a significant enhance- 
ment is observed for the contribution originating from 
the transition N — > 7rA which is not Pauli suppressed. 
These results emphasize the importance of careful treat- 
ment of the in-medium isobar self energy and propaga- 
tor, which is made possible by the convenient complete 
basis introduced in this context in Ref. [3. As a conse- 
quence in the nuclear medium the combined effects from 
the pion and from the isobar can produce sizeable in- 
crease in the pion light-cone distribution. The latter is 
constrained to smaller light-cone-momentum ratio y val- 
ues because of kinematical effect of the isobar-nucleon 
mass difference but can still have significant effects on 
the DY cross-section ratio. 

Since we are considering isospin symmetric nuclear 
medium and make a comparison with the deuteron it 
is advantageous to consider the pion distribution in an 
"isoscalar" nucleon, i.e. to consider a proton-neutron av- 
erage. Taking into account pions of all charges gives the 
complete pion distribution of an "isoscalar" nucleon: 



r /A (y) = 3/ 



TT°N/A 



{y) + V 



it A/A 



(y)- (36) 



In Fig. [5] we show the function f*/ A (y) for diferent input 
parameter values compared to the pion distribution of 
the free "isoscalar" nucleon. The probability Za of the 
bare nucleon in the medium takes the values from 0.6 to 
0.65, i.e. just slightly smaller than in the free nucleon 
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FIG. 7. (Color online) Isobar contribution to the pion dis- 
tribution (f A ^ A (y)) in the in-medium nucleon compared to 
the free nucleon. For both in-medium curves M* = 0.89 GeV 
and E V N = 0. 



FIG. 9. (Color online) The u p mcdinm) /u p bcc) (solid lines) 
and dp mcdmm '/dp frcc ' (dash lines) ratios of antiquark distribu- 
tions of an in-medium proton relative to the free proton for 
parameter sets (1) and (2). 



free nucleon 

M.=0.8GeV, g'„=1.0, g',^0.4, g' a =0.3 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

y 



FIG. 8. (Color online) Full pion distribution (f* /A (y)) in the 
in-medium nucleon compared to the case of free "isoscalar" 
nucleon. 

Before examining the DY cross-section we show the 
ratio of the antiquark distribution in the in-medium pro- 
ton and the same distribution in the free proton. The 
up and down antiquark distributions experience different 
in-medium modification due to different weights of nu- 
cleon and delta contributions even in isospin-symmetric 
nuclear medium. In Fig.|9]we show the ratios of antiquark 
distributions for an in-medium proton relative to the free 
one for two typical parameter sets denoted by "pion (1)" 
and "pion (2)", already used for plots in Figs. [3] and 
HI We observe pronounced enhancement for the up an- 
tiquark coming from the substantial pion enhancement 



due to the delta-baryon final state as compared to quite 
modest enhancement and even suppression for the down 
antiquark as a consequence of larger weight of nucleon fi- 
nal state and smaller weight of delta-baryon final state as 
compared to the up antiquark. This difference points to 
the possibility of distinguishing between effects coming 
from the medium modification of the nucleon and delta 
baryon by examining observables to which up and down 
antiquarks contribute with different weights. 

We now turn to the DY cross-section ratio (|10p . In 
Figs. [T0l and [TT1 we show the cross-section ratio (ITU1) as a 
function of x% for fixed values of X\ . The input parameter 
values are given in the figure caption. We observe an 
enhancement only for small values of X2, typically less 
than 0.2, and for X2 > 0.1 a decreasing trend as a result 
of the convolution with nucleon distribution (I17[) . 

For comparison with the measurements of Ref. @ we 
computed the ratio of the nuclear and deuteron cross 
sections for given Xi and integrating over x\ satisfying 
the condition x\ > X2 + 0.2 corresponding to the exper- 
imental cut-off. Fig. [T2] shows the measured values with 
error bars and the calculated curves for different input 
parameters. We consider the lowest curve in Fig. 1121 with 
M* = 0.8 GeV and corresponding rather pronounced sup- 
pression of the order (M*/M) 2 probably exaggerating 
the effect of the nucleon mean-field approximation and 
regard the other two curves as representing better our 
results based on the preferred parameter sets. 



VI. SUMMARY 

In this work we presented an analysis of nuclear effects 
on the Drell-Yan process. The approach is based on the 
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FIG. 10. (Color online) The cross-section ratio (|10[1 as a func- 
tion of x 2 for fixed values of xi. The used parameter values 
are: M, = 0.85 GeV, Ejv = 0.04 GeV, E^ = -0.1 GeV, E^ = 
0,^! =0.8, g[ 2 = 0.3,g' 22 =0.3. 



FIG. 12. (Color online) Experimental results from Ref. [jj 
compared to our calculation for different parameter values. 
Short dash line: M, = 0.89 GeV, E^ = 0; solid line: 
M, = 0.85 GeV, E^- = 0.05 GeV; dash-dot line: M* = 
0.8 GeV, Yl v N = 0.09 GeV. For all three curves: E^ = 
-0.1 GeV, E^ = 0, g' n = 1.0, g' 12 = 0.4, g' 22 = 0.3. 
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FIG. 11. (Color online) The same as Fig. \W\ but with 
used parameter values: M* = 0.89 GeV, E^ = 0, E^ = 
-0.1 GeV, E^ = 0, Pll = 1., g' 12 = 0.4, g' 22 = 0.3. 

pion-cloud model of the nucleon and a relativistically co- 
variant self-consistent in-medium calculation of the pion 
and delta baryon propagators taking into account nu- 
clear effects in the mean-field approximation. Starting 
with the free nucleon we showed that the observed d — u 
antiquark distribution can be well reproduced by suit- 
able choice of the 7rNN and 7rNA form factors with delta 
vacuum propagator taking into account its free width. 

Using the same values for the form factors we com- 
puted the pion light-cone-momentum distribution for nu- 
cleons in an isospin symmetric medium with density cor- 
responding to average densities of medium mass nuclei. 



We took into account the change of the pion cloud origi- 
nating from both the pion-nucleon and pion-delta states 
and a small correction (neglected in previous work) at- 
tributed to the binding effect of bare nucleon. 

Fermi motion and binding of nuclear nucleons were ac- 
counted for by the two-parameter light-cone-momentum 
distribution (TT71) which reproduces the negative slope of 
the classical EMC effect in the region 0.1 < x < 0.5 
as shown in Fig. [J] Taking into account the pion en- 
hancement which comes from the pion-delta state of the 
nucleon (which is significant only for small light-cone- 
momentum (y w 0.2) values) leads to some enhancement 
of the F2(x) ratio for x < 0.2. 

Pion and delta properties in the nuclear medium are 
calculated in a recently developed fully covariant self- 
consistent model [HI which consistently takes into ac- 
count the 7rNN and 7rNA vertex corrections due to 
Migdal short-range correlations. Pronounced softening 
of the in-medium pion spectrum present in simpler mod- 
els does not appear in this approach and consequently 
Pauli blocking causes some suppression of the pion dis- 
tribution coming from the pion-nucleon state for an in- 
medium nucleon. However, enhancement results from a 
careful treatment of the pion-delta state as a consequence 
of pion broadening and delta shift and broadening. 

The net effect for preferred parameter values is a mod- 
est enhancement of pion light-cone-momentum distribu- 
tion mostly concentrated around y rs 0.2 value. As a 
consequence the DY cross-section ratio exceeds one for 
small values, typically less than 0.2, of the X2 variable 
for fixed x\ values or integration over it corresponding to 
some experimental cuts. The convolution with distribu- 
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tion (fT7|) acts qualitatively on the antiquarks in the same 
way as on the quarks producing a negative slope which 
is less pronounced for the integrated cross section. 

As one can see in Fig. [T^] the large error bars of 
the measured cross-section ratio do not allow a sensi- 
tive enough comparison with calculated results in order 
to determine more precisely the preferred parameter val- 
ues and consequently the in-medium properties of the 
pion and delta baryon. It would be very desirable to 
achieve measurements with considerably smaller uncer- 
tainties which could contribute to the resolution of some 
decades old issues of nuclear physics. We remark that 
another interesting possibility for studying sea quark dis- 
tributions in nuclei would be the use of an Electron-Ion 
Collider as detailed in the joint report of the Brookhaven 



National Laboratory, the Institute for Nuclear Theory 
(Seattle, WA) and the Thomas Jefferson National Accel- 
erator Facility [3l|. The proposed semi-inclusive deep- 
inelastic electron-nucleus scattering would provide new 
information about the structure of nuclei and quantum 
chromodynamics of nuclear matter and extend possibili- 
ties for studying effects of the transverse momentum dis- 
tribution of partons [32j . 



Appendix 

The terms of expression (f33|) necessary to calculate the 
contribution of the in-medium delta are given as follows: 



with 



v4(A) 
99 

A w 

qu 

AM 

UU 

JQ) 



= 2 



-~qq 
~(Q) 

^qu 
-'qu 

,(Q) 



1 + A A + q-u B A /y/q 2 - (q ■ u) 2 

2 Re 



9 2 5a/vV- {q-u) 2 {l + A A -q-u B A /^q 2 - (q ■ uf) 

2 



q 2 B A /^q 2 -{q-uf 



(Ah + (p - T? N u) -P') t+(q- p'f - (q ■ X{p')f 



~(M t + (p - • p>) [t+(q- p'f + 3 (q ■ X(p')f 

-(M* + (p--E v N u)-p')(q-u -q-p'u-p' + q ■ X(p') u ■ X(p')), 
~\{M* + {p- Z v N u)-p'){q -u -q.p'u-p' - 3g • X(p')u- X(p')), 



ci Q J = (M. + (p - £» • p') [-1 + (u • p'f (u ■ X(p')f 



2 (M* + (P - E» ■ f>') [-1 + (u ■ P'f + 3(« • X(p')f] , 



a a ^ [nw 5 w (i - nWjW)- 1 ]^ + (i - n^ g ^)-% 3 , 

b a = [nw 3 w (i - nV)-^ + [n (i) 5 (i) (i - n^ g ^)-% 3 , 



x( P ) 



_ (P ■ u) p p - p 2 



P 2 y/{p ■ U) 2 /P 2 



(A.l) 



(A.2) 



where t = —q 2 and p M — p^/ * 



1 p 2 . 
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